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We investigate the exact solution, perturbation theory and master equation of open system dy- 
namics based on our serial studies on quantum mechanics in general quantum systems [An Min 
Wang, quant-ph/0611216, quant-ph/0611217]. In a system-environment separated representation, 
a general and explicit solution of open system dynamics is obtained, and it is an exact solution 
since it includes all order approximations of perturbation. In terms of the cut-off approximation of 
. . . perturbation and our improved scheme of perturbation theory, the improved form of the perturbed 

' solution of open systems absorbing the partial contributions from the high order even all order 

' approximations is deduced. Moreover, only under the factorizing initial condition, the exact master 

, equation including all order approximations is proposed. Correspondingly, the perturbed master 

equation and its improved form different from the existed master equation are given. In special, the 
^ ' Redfield master equation is derived out without using Born-Markov approximation. The solution of 

Q , open system dynamics in the Milburn model is also gained. As examples, Zurek model of two-state 

open system and its extension with two transverse fields are studied. 
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I. INTRODUCTION 



in 
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A realistic quantum system is never isolated, but is immersed in the surrounding environment (alias bath, reservoir) 
I and interacts continuously with it. Such a system without ignorable coupling to the environment can be called open 
\^ ■ quantum system. Generally, the environment consists of a huge number of degrees of freedom, it is even whole outside 
■ world (universe) of the concerning open quantum system. In fact, we might not know the exact state of the outside 
,: world, having only some statistical information to describe it. However, we are really interested in a reliable and 
effective theory of open system dynamics under the influence of its environment. 

The basic idea of quantum theory of open systems is thought of as an interesting open system and its surrounding 
(~| ' environment form a total composite system, or vis versa, a composite system can be decomposed into an interesting 
Ci ' open system and a surrounding environment. The key matters of quantum theory of open systems are to determine 
^ ' the interaction between the open system and its environment and build the physical models of the open system and 
~J ' its environment. Open system dynamics is just a law of this open quantum system how to evolute with time and its 
^ . solution at any given time. 

' Open quantum system and its dynamics are very important for many interesting quantum theory branches such as 
ijjj , quantum optics [l|, 0] , condensed matter theory, quantum information and computing Q , more concretely, quantum 
d ' decoherence, quantum measurement @,|6[, quantum dissipation 0, quantum transport @, quantum chaos et. 
al. Study of open system dynamics is helpful for understanding some very essential problems in physics, for example, 
the transition from quantum to classical world. 

A variety of different formal techniques have been developed and used in dealing with open quantum systems. 
From the above reviews and books, the interested readers may get them. Here, we intend to start with "the first 
principle" of quantum mechanics, that is, the Schrodinger equation or the von Neumann equation, and then try to 
build a theoretical formulism including the general and explicit forms of motion equation, dynamical solution, and 
perturbation theory of open systems. 

It is clear that such a "first principle" scheme might be not suitable to the cases when one cannot clearly know 
the environment model and/or the system-environment coupling form since the environment is too huge and too 
complicated. However, our conclusions might be helpful for building the models of such some systems. Moreover, one 
of possible ways to avoid this difficulty is to use the Milburn dynamics [ll|. That is, the environment is separated 
into near- and remote two parts, the Hamiltonians of the near environment (often with finite degree of freedoms) 
and the coupling to the interesting open system are assumed to be clearly known, but the influence of the remote 
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environment on the interesting system is incarnated by an extra term in Milburn motion equation compared with 
the von Neumann motion equation. Similarly, we successfully obtain the general and explicit solution of Milburn 
dynamics of the interesting system according our scheme. 

Because of the dissipative nature of open systems, we must turn to the density matrix for a proper description 
whatever the initial state is pure or mixed. Actually, we are interested in the properties of open systems only, it will 
be appropriate to study the reduced density matrix evolution with time or motion equation or its solution. Here, 
the reduced density matrix describing the open systems is obtained by tracing out the degrees of freedom of the 
environment from the total (system plus environment) density matrix. 

Due to the system and environment being entangled generally in system evolution with time, directly solving 
Schrodinger equation or von Neumann equation of the total system is a formidable task by using the existed methods. 
Traditionally, this problem is studied by perturbation theory in system-environment coupling scheme. Ones often take 
the interaction between the system and its environment as a perturbed part and then use the interaction picture to 
derive out the master equation of open systems via some physical approximations such as Born-Markov ones and the 
others. If an open system is exactly solvable, the coupling Hse is weak, the evolution time is short enough, and the 
used physical approximations are indeed appropriate, this has been proved to be an effective method. However, when 
the above conditions do not satisfied sufficiently, the problem gets complicated and perhaps leads to some difficulties, 
although some formal techniques have been developed and used in order to overcome some possible shortcomings. 
For generality and reliability in theory, we feel that we have to consider whether these approximations are necessary, 
if without these approximations, can we obtain the formulism of open system dynamics? The conclusions obtained 
here answer these problems. 

In this paper, we will provide the amelioration of the existed scheme of open system dynamics and try to build a 
theoretical formulism using our recent investigations on quantum mechanics in general quantum systems |12| . [Tsj . We 
first obtain the exact solution of open systems including all order approximations of perturbation and then give the 
improved form of perturbed solution of open systems absorbing the partial contributions from the high order even all 
order approximations of perturbation. Only under the factorizing initial condition, we derive out the exact master 
equation and its perturbed form via the standard cut-off approximation of perturbation. Moreover, we propose the 
improved form of perturbed master equation. In special, based on our master equation, we re-deduce the Redfield 
master equation without using Born-Markov approximation, and we point out the differences between our master 
equation and existed ones. We also obtain the solution of open system dynamics in the Milburn model. In order 
to illustrate our open system dynamics, we study Zurek model of two-state open system and its extension with two 
transverse fields. We are sure that our open system dynamics can be used to more open systems since its generality 
and clearness, and its calculations are simpler and more efficient, its results are more accurate and more reliable than 
the existed scheme. 

This paper is arranged as follows: besides Sec. [J is an introduction, in Sec. |TT1 by virtue of a system-environment 
separated representation, we first obtain a general and explicit solution of open systems including all order approx- 
imations; in Sec. Illli we gain the improved form of perturbed solution of open systems, which absorbs the partial 
contributions from the high order even all order approximations of perturbation; in Sec. IIVI we deduce the exact 
master equation of open systems only under the factorizing initial condition; in Sec. |V]we get the perturbed form of 
our master equation and its amelioration; in Sec. IVIl based on our master equation, we re-deduce Redfield master 
equation without using Born-Markov approximation, and we point out the differences between our master equation 
and existed ones; in Sec. IVIIl we obtain the solution of open system dynamics in the Milburn model for the Milburn- 
type closed total-systems. This implies that our solution and methods are applicable to more general open systems; 
in Sec. IVIIIl we study Zurek model of two-state open system and its extension with two transverse fields; In Sec. IIXI 
we summarize our conclusions and give some discussions. 



II. GENERAL AND EXPLICIT SOLUTION OF OPEN SYSTEM DYNAMICS 



In this section, we will derive out a general and explicit solution of open systems by using our recent work of exact 
solution in general quantum systems p^ . 

As is well-known, if assuming that the interesting open quantum system and its environment are taken as a closed 
(or isolated) larger composite system, that is, a total system, we can think that this total system obeys the Schrodinger 
equation or the von Neumann equation, respectively, for a pure state |4'fot(i)) or a mixed state ptot{t), that is 

-i^|^*ot(i)) - Htotl'ftotit)), (1) 

ptotit) = -i[Htot,Ptot{t)]. (2) 
where the total system Hamiltonian Htot that we consider here is made of the sum of the interesting open system 
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Hamiltonian Hg and its surrounding environment Hamiltonian plus an interaction i?sE between the system and 
the environment, that is 

Htot = Hs + He + HsE- (3) 

Note that the total system Hilbert space Htot is defined by the direct product T^s ® He of open system Hilbert space 
T^s and its environment Hilbert space He- Here and in the following, we will discuss time- independent Hamiltonian 
and we have taken h = 1 for simplicity. 

In an open system dynamics, a key difficulty to lead to the problem becomes intractable is that there is the 
interaction between the open system and its environment with huge degree of freedom. With time evolution, the 
open system inevitably entangles with its environment. Therefore, we starts from a system-environment separated 
representation (SESR). This representation is beneficial for obtaining the general and explicit solution of open system 
dynamics as well as proposing the improved scheme of perturbed theory [13i] , because in the SESR we can conveniently 
trace off the degree of freedom of environment. Introducing the SESR is a simple and natural idea, and we will see it 
is also very useful. To this purpose, we first divide the Htot into two parts 

Htot — Htoto + Htoti, (4) 

and, without loss of generality, we denote 

Htoto — Hso + Heo + HsEo, Htoti — Hsi + Hei + Hsei- (5) 

It is clear that 

Hso = hso <^ Ie, Heo /s ^eo, (6) 
while we need the coupling Hamiltonian with the following form 

^^SEO = ^'^CmnSmO ^ BnO- (7) 
m,n 

It is general enough if we do not restrict the forms of Smo and i?„o- In the above expressions, the total Hilbert space 
is Htot — Hs He, Is and Ie are, respectively, the identity operators in Hs and He, and Cmn are coupling constants 
between the open system and its environment. Note that -ffso and Heo are always hermitian as usual. In addition, 
we need Hseo be also necessarily hermitian. In fact, because Hso and Heo commute, the SESR always exists. The 
aim to add Hseo is to obtain better precision and to simplify the perturbed part when passing to perturbation theory. 
It must be pointed out that the general principle to divide H into two parts is to let the terms as more as possible 
belong to Hq but the precondition is that there exist the commuting relations: 

[hsO,Smo], [^EO, Cmrti^no] — 0, Or [hsQ, CmnSmo], [hEO, Bno] = 0- (8) 
n m 

Moreover, that the eigenvalue problem of Htoto is solvable. In fact, this solvability implies that hso and Heo are 
solvable, then hso and Sm, ^eo and J^n'^mnBn have the common eigenvectors, or hso and ^jy^CmnSm, hEo and 
have the common eigenvectors i.e, we have, respectively, 

hsolf) = E^\(l)'^), Smolf) = Sm-fl'l)''), ^EolXv) ^ £v\Xv) , ^ Cm7iBno\Xv) ^ rmv\Xv) (9) 

n 

and 

^S0|<^''> = E^lcf)'^), ^CrnnSrno\<l>'^) = Sn',\<l>'^), hEo\Xv) = £v\Xv) , Bno\Xv) = rnv\Xv) ■ (10) 
m 

They indicate that the eigenvectors of Htoto ^ or the common eigenvectors of -f^so, -f^Eo £^nd Hseo ^"^^ 

\<S>y-)^\^f)(g>\x-), (11) 
which span a separate representation of the system and the environment, and it is clear that 

Htoto\'i>^'')=E^,\<^^-), (12) 
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(13) 



It must be emphasized that the principle of Hamiltonian spUt is not just the best solvabiUty in more general 
cases. If the cut-off approximation of perturbation is necessary, it requires that the off-diagonal elements of the 
perturbing Hamiltonian Htoti matrix in the SESR is small enough compared with the diagonal elements of Htot = 
Htoto~^ Htoti matrix in the same representation according to our improved scheme of perturbation theory. In addition, 
if there are the degeneracies, the Hamiltonian split is also restricted by the condition that the degeneracies can 
be completely removed via the usual diagonalization procedure of the degenerate subspaces and our Hamiltonian 
redivision, or specially, if the remained degeneracies are allowed, it requires that the off-diagonal elements of the 
perturbing Hamiltonian matrix between any two degenerate levels are always vanishing, in order to let our improved 
scheme of perturbation theory work well I13|. As an example, it has been studied in Sec. IVIIII 

From the formal solution of the von Neumann equation of the total system 



Ptot{t) = e 

and our expression of the time evolution operator [l 



ptot{0)e 



iHtatt 



-iHtott 



1=0 

it immediately follows that the solution of total system density matrix with time evolution is 

00 00 
Ptotit) = Mt)Ptot{0)M-t) = J2 Ak{t)ptot{0)Aj{t). 

k,l=0 k,l=0 

In the SESR, we have 

Ptotit) = E E E Afit)p^-^p'-'mf-''''''^'{^twn(^'''^' 

I3,u,j3' .u' 7,v,7',d' k.l — 1 
00 

= E E E^fe''w^''"''''"'(o)<"'"''''(-t)[i0^ 

0,u,l3' ,u' 'y,v,'y' ,v' k,l — l 



where 

pZo^t^'^^) = ($^1Ptot(0) 

In Ref. [l2|, we have found the explicit forms of Ai{t). In the SESR, they read 

'it) 



(j)7 V 



^0 



(14) 
(15) 

(16) 

(17) 
(18) 

(19) 
(20) 

(21) 



E E 



'l+l 

E(-ir^ 

.i=l 



toil 



d,iE[-,v,l]) 

^7i7^'Ui'U^7i_|.i7'^'Ui_|_ii;' : 



(22) 



toil 



i^fi^i form so-called "perturbing Hamiltonian matrix", that is, the represen- 



and all iJ, 

tation matrix of the perturbing Hamiltonian in the unperturbed Hamiltonian representation (SESR). While 



(-^71"! -^7i + li'j=l) 

i=l 
i-1 

n (^7.'', 



diiE[jvJ]) 
d,iE[jv,l]) 
di+iiE['yv,l]) = n (^7,^, - -^7i+i^i+i 



E. 



} n (^7.^. 



(23) 
(24) 
(25) 
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By tracing off the degree of freedom of environment space, we obtain the exphcit expression of time evolution of 
reduced density matrix of open system 



Syyi , which is an advantage of the SESR. 



where we have used the fact TrE ^ = \il''^)(^ip^' 

It is clear that in the above expression, we need to know the concrete forms of {(f)'^) and \x^) in order to obtain the 
explicit expressions of Al^'^^(t). In fact, this is a physical reason why we take the form of Htoto as Eq.© so that the 
eigenvectors and eigenvalues of Htoto are obtainable. 

Note that there are apparent divergences in the above exact solution. For the tidiness in form, we keep these 
apparent divergences in our expressions, but we can completely eliminate them by the limit process [T§|. In other 
words, our exact solution of open systems should be understood in the limitation sense. 

Just as pointed out above, there is, at least, an inherent SESR (ISESR) in the total system if taking Htoto = Hs+H-^. 
We will be able to obtain the similar solution as Eq. ([26]). However, the ISESR is not unique in general because, in 
principle, a part of Hg and/or a part of He can be absorbed to Hse if -^s and /e are thought of as, respectively, the 
system operator and the environment operator. In this sense, the difference between the SESR and the ISESR is that 
the SESR allows to contain a part of Hsm = Sm.n <^mnSjnO B^o, in which, S'„o 7^ h and J^n Cmn^no 7^ -^e for all 
m. Of course, if the cut-off approximation of perturbation is necessary, it requires that the absorbed parts from Hg 
and He arc small enough. Such an example is discussed in Sec. IVIIII In addition, one of the reasons to introduce 
the SESR is to make the Hamiltonian redivision and absorbing the perturbing parts of Hg and He to the perturbing 
Hamiltonian of the total system look like more natural. 

Different from the general and explicit solution HH), the coefficients of our above solution (|26p of open system 
dynamics are c- number functions whose forms are expressed clearly. Because A2^'^^{t) include all of order approx- 
imations, this solution is, in fact, exact although it is an infinite series. Our solution in form is general enough, 
and it is able to applied to the cases that Hg and/or He are not exactly solvable. It is clear that we do not use the 
accustomed approximations such as the Born-Markov approximation, the factorization assumption for the initial state 
et. al. Hence, it should be more general and more reliable in theory. Moreover, by virtue of the improved scheme of 
perturbation theory proposed by us, we can obtain the improved perturbed solution of open system dynamics with 
better precision and higher efficiency because the contributions from the high order even all order approximations of 
perturbation are absorbed into the lower order approximations. 



III. IMPROVED PERTURBED SOLUTION OF OPEN SYSTEM DYNAMICS 



Traditional scheme of perturbation theory has been successfully used to solve many systems. However, in our point 
of view, it is still improvable, even it has a flaw because it introduce the perturbing parameter too early so that 
the contributions from the high order even all order approximations of the diagonal and off-diagonal elements of the 
perturbing Hamiltonian matrix are, respectively, inappropriately dropped and prematurely cut off. For some systems, 
the inffuences on the calculational precision because of this flaw can be not neglectable with the evolution time 
increasing. Actually, the traditional scheme of perturbation theory does not give a general term form of expanding 
coefficient evolution with time for any order approximation and does not explicitly express the general term as an 
obvious c-number function. Thus, it is necessary to find the perturbed solution (or perturbed energy and perturbed 
state vector) from the low to the high order approximation step by step up to some order approximation for a needed 
precision. Recently, in our work, we pro pos ed an improved scheme of perturbation theory based on the general and 
explicit form of our exact solution |12l . [l3|. In our improved scheme, we introduce the approximation as late as 
possible, and consider subtly and systemically the affection of high-order approximation to the low-order one by the 
dynamical rearrangement and summation method. This flnally results in the improved form of perturbed solution, 
and its expansion coefficients contain reasonably the high-order energy improvement. In this section, we will apply 
our improved scheme of perturbation theory to open systems. 

It must be emphasized that before applying our improved form of perturbed solution, we have to first carried out 
the digonalization of degenerate subspaces if there is degeneracy and do the Hamiltonian redivision when Htoti bas 
the diagonal elements, in order to completely removed possible degeneracies by this procedure. When the remained 
degeneracies are allowed, it requires that the off-diagonal elements of the perturbing Hamiltonian matrix between any 
two degenerate levels are always vanishing. For more complicated cases, we will study in the near further. 
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Therefore, up to the three order improved approximation, we have 

3 



fc+i<3 



(27) 



where 



e ^ 



E^y E-y'yf E^y E-yfy' 



9{ 



(28) 
(29) 



A]^'^ " {t) 



I „7«'.7i''i „7ii'i,7''A ,A , -I- 

71,1)1 L yt^iv -t^J-Yivi) 



lE^lylt 



e T 



i^E^y -^'yiUl) (-^71U1 -^7'lj') (-^7U Ery'y/^ (^E^^y^ E^'yl^ 



7tJ,7ii;i 7it)i,7 D 1 

i/i yi il-iv,j'v' ( 1 



(30) 



^7^.,7'.'(^) = E 



i/l i/l Ul U-y^^'Oyy' 



{E^y E^^y^ ) i^E^^y^ "^72 ^^2 ) (-^7V ■^J2V2) (^71^1 ^llV2) 



-E 

71 

+ E 



i^E^y E^^y^) (^E^y Ery'y'^ (^E^y E^^y^^ i^E^y E^'y'^ 



7i>72 L 



-'77; -'-'71 1)1/ V 7^^ j-^^'' 



7t;,7ii)i 7ii)i,7t; 7^7 

i/i 5i ffi 



e ' ^1^1 ij-iivi.-f'' 



i^Ejy Ej-^y-^^) {E^y E^^y^^ {E^y E^'y') {Ejy E^^y^^ (^E^^y^ E^^y^^ (^E^^y^ E^'y'^ 



Q ^^l^^^^rj^y^^^y^ 



{Ejy -E^72'1'2 ) (-^7l'l'l -^72'"2 ) (-^72l'2 Ejlyl ) 



7t>,7iwi 71 j;i,72W2 72^2,7 " „ 
Ul ill ill iljv,j'v'j 



(31) 



where S^^' and Syy' are the usual discrete delta functions, while rj-y^i = 1 — ^-y^', riyyi = 1 — 6yy', and rj^y^j'y' = 
rj^^, + S^j'Tjyy' = rjyy' + Tj^j'Syy'. Moreovcr, we have defined so-called improved form of perturbed energy by 



E^y E^y G^y G^y G^y Giyy G^ y 



(32) 



where, G^^j* = arc diagonal elements of Htoti and g^'""^^''^ arc off-diagonal elements of Hton in the representation 
of Htoto- In addition, hj" include the diagonal elements after the diagonalization of degenerate subspaces. While 



= 1 



71 1"! 



Eryy Ery^ y^ 



7t>,7il)i „7lWl>7'' 

9i 91 ) 



(33) 



G'^^= E 



ll,V\,^2,V2 



{Ejy E^^y^^{E^y Ej^y^^ 



jv,-fivi 7i?;i,72t;2 72?^2,7V 

9i 9i 9i ) 



(34) 



^7V 



= E E 



71,72,73 Vl,V2,V3 



7ll,7lt)l 7lt)l,72f2 „72'"2,73''3 „73f3,7''„ 

9i 9i yi 9i '/7ii,72f2 

{Eryy E J -^y -^^ {^E^ y E^^y^^{E^y E^^y^^ 



EE 

71 >72 Vl ,V2 



7t>,7iJ)i 7ifi,7f 7i),72f2 „72'!'2,7'' 

9i 9i 9i 9i 

{Eryy Ey-^y^) {Eyy Ey^y^^ 



(35) 
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E 



E 



71,72,73,74 Ul,1)2,f3,l'4 



E E 

71 ,72 ,73 ''1 ,1'2 ,f3 



JV,JiVi 71 Dl ,72t)2 „72f 2 ,73^3 „73l'3 ,74l'4 74^4 „ „ 

yi yi yi yi yi 'n" ,72^2 '/7'",73"3 

7i>, 71111 711)1,711 7ti,72i'2 „72i'2,73i'3 _73i'3,7i' 

.'ii yi yi yi yi 

(^E^y E^-^y^^ (^Ejy ^^2^2^ {-^fV "^73 '1'3 ) 



7«,7ii,i 7iDi,7i' 7i),72f2 „72'"2,73'J3 „73i'3,7'' 

yi yi yi yi yi 



7'U,7it)i 7if i,7f 7i),72f2 „72'u2,73i>3 73i'3,7'' 

i/i i/i yi i/i yi 



{E^y Ej-^^y^^ (^E^y Ej^y^') i^E^y E^.^y^^ (^Ejy E^-^y-^^ (^Ejy E~y^y^^ (^Ejy Ej^y^^ 



(36) 



It must be emphasized that if only based on our calculations that was completed in Ref . |13l | , the improved perturbed 

energy in the exponential powers of An, Au and A13 will be cut off, respectively, to G^Ju-, G^Jy- and g'^Ju-. However, 
according to our conjecture, we think that they can congruously written as the definition ([5^ . 

Our improved perturbed solution inherits some features from our exact solutions, for example, it is an explicit 
c-number function, easy to calculate, docs not need the extra approximations. In principle, we should can calculate to 
any order of improved approximation. It must be emphasized that our improved form of perturbed solution absorbs 
the partial contributions from the high order even all order approximations of perturbation. This means that our 
solution has better precision and higher efficiency. In fact, these advantages have been seen in our recent work p!^.[l3|. 

IV. MASTER EQUATION OF OPEN SYSTEMS INCLUDING ALL ORDER APPROXIMATIONS 



Because we have obtained the general and explicit solution of the open system dynamics when the Hamiltonians of 
the system, its environment and the interaction between them are known, it is unnecessary to derive out the dynamical 
equation of open systems. However, in order to understand the affection from the environment, compare our solution 
with the existed motion equations and reveal the improvement of our method, we would like to discuss the motion 
equation and master equation in this section. 

It is more convenient to derive out the master equation in the inherent SESR (ISESR) of open systems. That is, 
we take Htot = Hg + H^. In fact, it make us more easily compare our results with the existed ones. Obviously, the 



bases of ISESR are IV''' 



that is 



Hs\ib"')(g>\Lo-)^Es^\^^)(g>\to^), 



(37) 
(38) 



Similar to the way in Sec. |IT1 we can obtain the exact solutions ptot{t) and ps{t)- AH we need to do is to change 



r-i' : ' ' ' 

Ix") as ® and define the all matrix elements in the ISESR, for example, yl^ " 'T ^ = {^"^ uj'"\Ahif) 
Therefore, 



Ptot{t) 



E 

A;,^— I3,u,(3' ,u' 7,1), 7',?; 



(39) 



^sW = EE E Ar''"Wefo"/"'(o)<"''''"(-t)iV'-^)(^^' 

7?;, 7' k,l—0l3.uJ3'.u' 



(40) 



From the solution ([39|) . it is easy to get that TrE {[Hs , ptot{t)]} = [/is, PS (01 ^^i^ TrE {[He, ptot{t)]} = 0. Hence, 

ps{t) = TYEPtotit) = -iTrE {[Htot, Ptotit)]} = -i [hs,Ps{t)] - iTrE {[Hse, Ptot{t)]} , (41) 
where /is — Tt^eHs- Denoting system operators by Sm and bath operators by B'^, the most general form of Hse is 

Htotl =HsE^Yl ernes'™ (g) = ^Sm® Bm, (42) 

m,n m 

where Bm = X^n*^™"^"- Substituting the above relation into Eq. (|4ip . we obtain the motion equation of open 
systems 



ps{t) = -i[/is,Ps(i)] -i^[^m,TrE{(/s®B„)ptot(0}]. 



(43) 
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The second term of its right side represents the influence of the environment on the system. 

In order to express the motion equation l|43p of the open systems in the exphcit matrix form, we introduce so-caUed 
factorizing initial state assumption, that is, the system and its environment are uncorrelated initially such that the 
total density matrix is a direct product of the system and its environment density matrices. 



Ptot(O) =ps(0)®pe(0). 



(44) 



Its advantage is to make us easily consider the actions of the operators on, respectively, the open system space and 
its environment space, and finally we can easily trace off the degree of freedom of environment space. In order to use 
this advantage, we introduce two new operators 



(45) 
(46) 



1,1' 



where Vs{P, /?') = (jjj^ are the basis operators of the system Hilbert space Hs, while the operators A^}^{t, j3, fi') 



and ^g'^2(— t, 7, 7') are defined in environment Hilbert space TYe by 



u.u' 



(47) 
(48) 



Thus, we see that A^l^'^ (t) and {—t) are decomposed as the summations whose every terms with the form that the 
open system parts and its environment parts are separate. Hence, we obtain 



Ptot{t)= J2 [^s(/3,/3')^?s(^)^s(7,7')]® [-4|3'l(t,/3,/9')eE(i)^l3'k-^:7,7') 



where 



0s{t) 



ps(0)e-"'^*. 



QEit) - e-»'-VE(0)e-'''-*, 

and then 

gtotit) = Aoit)ptot{0)Aoi-t) = gsit) (E) QEit). 
Substituting Eq. into the motion equation it immediately follows that 



(49) 



(50) 
(51) 



(52) 



ps(i) = -iN,Ps(t)]-iE E E ^r/^"yW[^"'^s(/3,/3')fs(<)P(7,7') 

m fe,i=0/3,/3',7,7' 



where we have used the fact that 



C^^!^^,{t) = TrE \B„,A'il{t,l3,(3')gE{t)A^^U^t,^,y) 



(53) 



(54) 



Further deduction needs us to obtain the motion equation of gs{t) that is expressed by ps{t)- In fact, based on Eq. 
(131]), we have 



gs 



it)^ps{t)- E E K','f,',,y{t)[VsiP,P')gs{t)Vsh,l')] 



k,l=0 P,0',l,l' 
k+l>0 



where we define the coefficients 



K^J.,^y{t) = TrE Uilit,f3,P')gEit)A^i^^i~t,^,j') 



(55) 



(56) 
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Therefore, we can use the iterative method to rewrite it as 



it) = 



M=l 



M 



"1=1 fc„,;„=o 0m.l3'^,7mn'„^ 

km+lm>0 



.1=1 

Substituting it into Eq. ([53|l . we obtain 



M 



M 

n^s(7,,7-) 



ps(t) = -i[/is,PsW]-iE E E C'r/.77'W['^™'^s(/3,/3')Ps(t)P(7,7')] 

oo oo A'' 

-iEE E ^"/;WwE(-i)" n E E 

m fc,i=0 /3,/3',7,7' 



Af=l 



AT 



N 



Sm,Vsi(3,(3') ( n^s(A,/3:) 1 Ps{t) I n^s(7j,7-) | ^s(7,7') 



(57) 



(58) 



Up to now, we have not introduced any approximation except for the factorization assumption for the initial state. 
Since our master equation l|551) including all order approximations, we can say it is an exact master equation of open 
systems. 



V. PERTURBED MASTER EQUATION OF OPEN SYSTEMS 

In the most cases, the interaction between the open system and its environment is weak. We can cut off the above 
exact master equation to some given order approximation. It is clear that since we absorb the coupling coefficients 



into Bm, we known C^'J^^^^, (i) is a quantity of the {n + I + l)th order approximation, KJ^^, ^^,{t) is a quantity of 
the (n + ^)th order approximation from their definitions. Although we can obtain any given order approximation of 
master equation based on our exact master equation (j58p . in most cases, we only are interested in up to the second 
order approximation. Because 



(59) 



^m,0/ r rp 

<-'/3/3',77' 



Trp 



A^l{t,P,P')QE{t)B^ 



(60) 
(61) 



we have 



where 



(62) 
(63) 



ps{t) = -i N,Ps(t)] - m),Ps{t)] ~iY.^S.r.,Psit)C^Rl{t) + ciZit)Ps{t)] 

in 

J(i) =^5,„(t)TrE [B^r^Q^it)), 



(64) 
(65) 
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Cill(i) = Tte { [h ® Bm] Ai (t)e'^--* [h QK{t)] } , (66) 

CnLit) = TVE{[/s<»eE(i)]e-^^-"UiM)[JsOB„]}, (67) 

LW(t) = TrE{yli(t)e'^*-°*[/s®eE(t)]}, (68) 

i?W(i) = T^E{[/s®^E(t)]e-^^-°Ui(-t)}, (69) 



while 



= ^ ^f"'^''(i)|^'5)(V'T|«)K)(a;''|. (71) 

Wc can sec that the Rcdfield master equation will be obtained from our this master equation without using Born- 
Markov approximation in next section. 

In order to absorbing the partial contributions from the high order even all order approximations into the lower 
order approximations, we can use our improved scheme of perturbation theory. In similar way used above, we have 

1 

ps{t) = -i [hs,ps{t)] + i [J{t), ps{t)] - i E E I-^-' Psi*)CiRm(t) + C[ll{t)ps{t)_ 

0=0 m 

-i [j{t),psm^\t)+Lf\t)Ps{t)]+\Y.[SrnAZ{t)Psmf\t) + C\^^^^^^ 

m 

+PBm^\t)diLit) + ^i'^(i)Ps(t)<L(t)] + 0{Hl,,), (72) 



where we have defined 



while 



CilL(*) = T^^{[h(^Bm]Aik{t)e'"'°'°'[h®m{t)]], (73) 

ctLii) = TrE{[/s®eE(t)]e-'^-''Un(-t)[/s®B™]}, (74) 

L':t\t) = TrE{^ifc(t)e'^--*[-fs0^E(i)]}, (75) 

Bil\t) = TVE{[/s®^E(t)]e-'^-°Un(-t)}, (76) 



^ii(*) = E '"'fT' _ '7^^^* gton^"'^" (1 - 50^5uv) ^ (77) 

= ^ ^f,"'^''(i)|V''>(^''|®K)(a;''|. (78) 



Here, E^y = Eg^ + s^y + h^y + G^y + g'^} + G^t) + • • • j E'^y = E^ + Sy + h^y, h^y are diagonal elements of i/, 



and the perturbed part of Hamiltonian in has be redivided as H^^^^ = Hse — ^jy ^7«|V'^)(V'^I <8i \aj"){ijj^\, that 

It must be emphasized that the operators in the above definitions and expressions are defined in the ISESR (that 
has been diagoalized in the degenerate subspaces if the degeneracy cases exist). However, Aik{±t) including E have to 
be calculated using H[ that is the perturbing Hamiltonian via the redivision skill. Hence, it is important to distinguish 
Htoto = Hs + -^E, Htoti = -ffsE and their redivision H[g^Q, fffoti in spite of them in the same ISESR. In addition, 
we assume all degeneracies are completely removed by the diagonalization procedure of degenerate subspaces and/or 
hamiltonian redivision for simplicity and determination. If the remained degeneracies are allowed, it requires that the 
off-diagonal elements of the perturbing Hamiltonian matrix between any two degenerate levels are always vanishing, 
in order to let our improved scheme of perturbation theory work well. 

In the above derivation of our master equation, we do not use Born-Markov approximation, but only standard 
cut-ofi: approximation. From our point of view, it is more reasonable in physics theory and its precision and reliability 
should be better in practical applications. 
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VI. RE-DEDUCTION OF REDFIELD MASTER EQUATION 



In order to compare our master equation with the known master equations and illustrate the validness of our 
master equation, we will deduce the Redfield master equation from our master equation without using the Born- 
Markov approximation in this section. In addition, we point out what differences between our master equation and 
the existed one, and provide the comments on well-known approximations using in open system dynamics. 

Firstly, we assume a thermal equilibrium for the environment, that is. 



Pe(0) 



-/Sb^e 



1 



(79) 



where /3b — ^/k^T with T the bath equilibrium. This is justified when the environment is "very large". Thus, it is 
easy to get 



^ ^ p-KEsf,-Es-y)t _ 1 ^ 



drTrp 







ihst 



where 



Likewise, we have 



-1 / drTrE [HsE{T)iIs(^ PEiO))] , 
/o 



(/s ® Hse{t) {Is ® ^^(0)) 

-i / drTrE [{Is <E> B„,) Hse{t) {Is Pe{0))] , 

^0 



(80) 

(81) 
(82) 

(83) 

(84) 
(85) 



drTrp 



{Is ^ Pe{0)) Hse{t) (/s®i?m(i) 



Jhst 



= i / drTrE [{Is Brn) Hse{t) {Is ® Pe{0))] ■ 
Jo 

Therefore, our master equation (|64p up to the second order approximation can be rewritten as 

ps{t) = ~i[hs,ps{t)]-i[J{t),ps{t)]- f drTrE { [i?sE, [i?SE(r),ps(0'»PE(0)]]} 

Jo 

+ J{t),J^ drTrE {[i?sE(T), PS (i)®PE(0)]] 



(86) 
(87) 



If we introduce the interaction picture, that is, an operator O in this picture is defined by a corresponding operator 
in the Schrodinger picture 



d{t) = e'^'°*°*Oe"'^'°*o*. 



(89) 



It is clear that for an operator Fs in the open system Hilbert space and an operator in the environment Hilbert 
space, we have 



Fs{t) = e"'s*i^se""'"*, 
FE{t) = e"''=*i^Ee-"'^*. 



(90) 
(91) 
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It immediately follows the master equation in the interaction picture: 



dps ft) 
dt 



-1 



J{t),Psit) 



drTr 







E { [-ffsE , [-ffsE (r) , PS ) PE (0)] ] } 



J{t), f dTTrE{[i/sE(T),psW®PE(0) 



(92) 



It must be emphasized that ps{t) is equal to e'''^*ps(i)e but not e'''^*ps(0)e 
In special, when we introduce the assumption 



Tr, 



{[i?SE,Ptot(0)]} -0. 



we have 



~"'«*TrE ■ 



Thus, Eq. ([55)1 becomes 



{[i?SE,Ptot(0)]}e"'^* = 0. 



(93) 



(94) 



ps(i) = -iN,Ps(0]-iE E E ^;^^WW[^™'^s(/3,/3')^>s(^)^(7,7') 



(95) 



™ k.l=0 /3,/3',7,7' 
fc+i>0 



that is, the perturbed form of master equation up to the second order approximation reads 

hit) = -iN,Ps] -iE[^™'/^s(i)C^lj'i(t) + C« (t)ps(t)]. 



(96) 



This means that the approximation ([93|1 leads to the following terms 



i[J(i),Ps(i)H 

or, equivalently, in the interaction picture 

-i[j(i),ps(<) 



J(i), / dTTrE{[i/sE(T),ps(t)'^PE(0)]] L 



J{t), 



drTr, 



{[i7sE(T),Ps(t)'»PE(0) 



(97) 



(98) 



are dropped by comparing Eq. with Eq. ([55]) or 

Usually, the approximation is thought of as a unimportant restriction since one can absorb the the dropped terms 
into the system Hamiltonian iJg. However, based on the above result, we think that the approximation ([M]) is a real 
assumption because the second term in (|97p or (|98p is not nontrivial and it can not be absorbed into i7g in general. In 
other words, the the second order contribution to the master equation from the second term in ([57]) or ([55]) should be 
considered and the approximation ([93]) should be rechecked for the concrete open systems except for the cases when 
J(t) = 0. Actually, we think, if J{t) is not equal to zero, the last term appears in our master equation or ([^^ is 
obviously different from the existed master equations. 

It is very interesting, when the approximation (j93p can be used to some given open systems, we immediately from 
the equation ((96)) obtain 



dps(0 
dt 



drTrT 



{[//SE, [i?SE(T),Ps(i)®PE(0)]]} 



(99) 



This is just the well-known Redfield master equation. This conclusion implies that the Rcdfield master equation is 
still valid without introducing Born-Markov approximation. Therefore, we think that Born-Markov approximation 
is unnecessary for the master equation with the second order perturbed approximations. From our point of view, 
this is a real physical reason why ones should use jointly Born- and Markov approximations and why ones can 
obtain useful conclusions in the cases without Born-Markov approximation. In fact, those terms that are dropped by 
Born approximation are compensated by Markov approximation. In other words. Born approximation plus Markov 
approximation back to no approximation based on our results. 



13 



VII. MILBURN DYNAMICS FOR OPEN SYSTEMS 



Historically, a useful dynamical model of open system is the Milburn model [ll| . It provides a way to describe so- 
called "intrinsic" decoherence. However, in our point of view, perhaps, it can be called external-external environment 
decoherence. That is, Milburn dynamics might be alternatively explained as the effect of environment of the composite 
system or the large environment of the proper system. This explanation is, in fact, a conclusion from that we believe 
the von Neumann equation is uniquely correct for a closed system. One argues what mechanism results in that the 
external influence is reflected by the extra term in Milburn dynamics. We can not answer it at present, but we 
would like to ask what condition changes the dynamics from the von Neumann's to the Milburn's. If the answer is 
the Milburn dynamics is a nature of closed quantum systems, then it is very difficult how to understand the free 
parameter ^o- 

Actually, only one can do something within the near environment in order to control decoherence, for example, the 
self- interaction of environment, and it is possible that one only knows how to appropriately describe the dynamics of 
near environment and the interaction between the system and near environment but are short of the knowledge about 
the remote environment. Therefore, in this section, we intend to use the Milburn model to consider the dynamics of 
the composite system made up of the system and its near environment. The conclusions obtained here imply the our 
solution and methods are also applicable to more general open systems such as the Milburn model. 

Dynamics in the Milburn model replaces the usual von-Neumann equation of the density matrix by 



Ptot{t) 



-\[Htot, ptot{t)] - —[Htot, [Htot, Ptot{t)]], 



(100) 



where is a constant meaning that there is some minimum unitary-phase transformation. This implies that coherence 
is destroyed as the physical properties of the system approach a macroscopic level. Hence, seemingly, the "intrinsic" 
decoherence explanation looks like to be reasonable. However, the minimum unitary-phase transformation is not 
clear. The parameter Oq in the Milburn model is still "free". In other words, 9q is not been given by the theory. If we 
think the extra term is resulted in by the remote environment, should be able be known by the experiment. 

Now, we directly extend Milburn dynamics to a Milburn-type closed quantum system consisting of the concerned 
system and its near environment. The Hamiltonian in eq. l|100|) still reads Htot = Hg + + Hse- Here, a Milburn- 
type closed quantum system is not really closed system from the view that a really closed system must obey the von 
Neumann equation. Actually, an alternative explanation is that a Milburn-type closed quantum system is still affected 
by the remote (larger environment), and this influence is represented by an extra term with multiplier because one 
cannot know the Hamiltonian of its remote environment and the interaction form between the interesting system and 
its remote environment. Obviously, when = 0, Milburn dynamics back to von Neumann dynamics. This implies 
that the (very) remote environment can be ignored. 

The formal solution of Milburn dynamics for the composite system can be written as (l^ 



Ptot 



{t) = exp{-iHtott - 0oHl,t/2} [e^Vtot(O)] exp {iHtott ~ OoHl^t/2} = Mkit)ptotiO)Mlit), (101) 



where 9Jl is a superoperator, i.e, dJlptot = (^oHtotPtotHtot, and the Kraus operators Mk{t) is in the form 



Mk{t) 



exp {-iHtott - eoHltt/2} 



(102) 



Without loss of generality, using of the denotation {A + B)^ — + /^(^, B) and f^{A, B) — 0, we can write 

oo n 

Kto exp {-iHtotot - eoHl,,t/2} 



Mk{t) 



\ r~im ^n-\-m-{-k ( TT j i \ 

] J [-titotO:-^totl) 



n—0 m— 



n! V 2i 



(103) 



Just as we find the exact solution of open system in von Neumann dynamics, we need a system-environment 
separated representation (SESR), which has been given in Sec. |TT1 Thus, in this SESR, based on the our expansion 
formula of operator binomials power [l3| we have 



K 



f''iHtoto,Hton) = E E E C;l'{E[^v,l] 

1=1 7l.--- ,71 + 1 vi,--- ,1)1 + 1 



n rr 7jf3,7j + l,i'j + l 
^totl 



(104) 
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i+i 



E. 



K 



d,{E[^v,l])' 



(105) 



where di{E['yv,l]) is defined in Sec. [TTl Therefore, the expression of Mk{t) is changed to a summation according to 
the order (or power) of the Htoti as follows 



M^,^'-' it) 



-^y — Eyyg{Eyu,t)d^^^>Syy' + 



kl 



E E E 

1=1 7i>--- ,7i + i 'Jir-- 



h d,{E[^v,l]) 



I 



^ J_ J_ -"toil 07i707i + i7'0'unj0t)i + it)' J 

J = l 



where the time evolution function g(x] t) with the exponential form is defined by 

g{x] t) — exp {— ia;t — 9i^x^t/2^ . 



(106) 



(107) 



Obviously, M^T"''^'"' (i) can be given via replacing g{x;t) by g*{x\t). Furthermore, we obtain the expression of time 
evolution of reduced density matrix of open systems, that is a general and explicit solution of open systems in Milburn 
dynamics: 

Psit) = J2 J2 Mf"'''"'(0/'"''^'"'(0)M2^'^'''^^<5„,|V'')(V'^| 



J2p^'"'''^iO)g{Ep.-E,y;t)\r){r\+Yl E /"'^'"'wE E 

/3,7,^' l3,u 'y' ,v' 7i,--- ,7^ + 1 



^ ' dk{E[jv,l]) 



.i=l 



E E 



^ ' dkiEWuJ]) 
1=1 L»=i ^' ' J'' 



P.P' ,u' 1,v 



00 00 



E E /'"'^'^'(o)EE E E 

P,u,/3',u' 7,t),7',t)' fc=l (=1 /3i,--- ,/3fc + i 7l.--- ,7! + 



\tp''){r\ 



k+1 l + l 

EE(-i)'"^' 

1=1 J=l 



^totl 



.i=l 



n rr 73f3,7j + l'"j + l 



d,iE[^u])d,iE[^v,l]) 



It is clear that if 6q = 0, this solution is just the form of solution of van Neumann dynamics that is obtained in Sec. 
ITTl Usually, the finite (even often low) order approximation about Htoti can be taken, thus this expression will be cut 
off to the finite terms. Similar to the methods used in Sees. IIIIl IIVI and IVl we can study the perturbed solution and 
motion equation of open systems in Milburn dynamics. It is not difhcult, so we omit them in order to save space. 



VIII. EXAMPLE AND APPLICATION 



In order to concretely illustrate our general and explicit solution of open system dynamics, we recall an exactly 
solvable two-state open system for decoherence that was first introduced by Zurek [l^ . In this Zurek model, the 
"free" (unperturbed) Hamiltonian TJso and the self- interaction (perturbing) Hgi of concerning two-state system and 
the 'free" (unperturbed) Hamiltonian H^q and the self-interaction (perturbing) H^i of the environment are taken 
as to be equal to zero. The total Hamiltonian of the composite system made of the interesting system plus the 
environment only has their interaction term, that is 



(109) 
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where the environment operator B^^ is defined by 

Ne /k-1 

= V 6^ /eJ ® {Zka'") (g, \ (y) /eH , (110) 





and A'e is the degree of freedom of the environment, which is very large even infinite. 

It is clear that this Zurek model can be exactly solved out. Its eigenvectors are so-called natural bases 



where ns,ni,n2, ■ ■ ■ = 0,1 and 



Insns) = |ns; ni, n2 • • •) = |ns) «) \nk), (111) 

k=l 

|0)s = (I), Il>s=(?)' (112) 

lA fo 



|o)fc = ( ' li^'^^ 1 )■ (113) 



The corresponding eigenvalues are 



EnsUE = Ensmn,- = (-1)"^ ^(-l)"'=Zfe. (114) 

fe=l 

Note that we use a simple notation to denote ni, n2, • • ■ here and after. 

Now, let we solve this Zurek model by using our exact solution or the improved form of the perturbed solution. 
From our point of view, the assumption that Hg and are taken as zero is a theoretical simplification. In fact, we 
can think that iJg and are constants so that we can absorb them into energy eigenvalues or, equivalently, directly 
omit them since these constants do not affect physics. Therefore, the base of the SESR can be taken as the natural 
bases (lllip . 

Since Hqe is completely diagonal in this SESR, that is 

(msmEl-ffsElnsnE) = (-1)"^ 2(-l)"^ZfeJ„3„3 [] (5™,„.. (115) 
We should use the Hamiltonian redivision skill, and then 

HtotO = Hzurek- (116) 

It is easy to get 

En.nE - -B„srun2... = (-l)"^^(-l)"'=^fc. (117) 

fc=l 

Au{t) - 0, (?>0). (118) 

where we have used the fact (^™smE,"snE _ q i^g^ggjj -ffjoji — 0. This means that the perturbed solution part of 
higher than the zeroth order approximation is vanishing. Therefore, our exact solution or the improved form of the 
perturbed solution including only non-vanishing zeroth order part becomes 

1 _ 

Pzurek{t) ^ e-'^"s."E*p™s™E,«s"E(o)giiJ„s"E*|^s^E)(nsnE| =e-'^-'--V(0)e'^-"'--*. (119) 

™Si"S=0 ITlE,"E 

Obviously, it is equal to the exact solution of the Zurek model (|109p via directly solving it. Of course, the solutions 
Ps{t) of this open system obtained by our exact solution or improved form of perturbed solution formula or directly 
solving method are consistent. Therefore, we can say our improved form of perturbed solution indeed absorbs the 
contributions from all order approximations of the perturbing Hamiltonian -ffsE since it is diagonal. In addition, we 
would like to point out that although there are the degeneracies in En^uE when mg + rrife ~ ns + nk, our improved 
form of perturbed solution can work well since H^^^i = 0. 
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In order to reveal the advantages of our exact solution and perturbed solution, we add two transverse fields, 
respectively, in the system and the environment, that is 

Htot = Mf^g ® /e + i^Zurok + ® S^E = AiCTg ® /e + (T^ ® (5:.E + S,e) , (120) 

where 



BxE — 



^ (g) /eJ ® [Xk^E") ® 
k \i=l ) \j=fc+l 



(121) 



The problem only with the system transverse field was studied in Ref. [IGj. The model (|120p is not exactly solvable 
unless /X = 0. Obviously, there are four kinds of the SESRs. 
Case one: The Hamiltonian split is 

HtotQ = Hs + Hb ^ fi(Js <S) Ie, Htoti^<Js^{B^E + B,E)- (122) 
The bases of unperturbed SESR are 

Ne 

IV's'x"^) = IV's')^(8)lx"'=), (123) 

where 

\r') = -L[|0)s + (-l)"^|l)s] (124) 
|X"^> = , ^ [{Z, + {-ir^Yu)\0),+X,\l),]. (125) 



where Yk = i/^f+Zj. Here, | X"*") are the eigenvectors of the environment operator B^e + B^e — (Xk<T^ + Z^cr^), 
and corresponding eigenvalues are (— 1)"*= Yfe. Thus, the eigenvalues of -ffjoto acting on IV's^X"'^) are 

Bjique — -^ns,nin2'-' — /^( f) ^- (126) 

Case two: The Hamiltonian split is the same as (|122p . and the corresponding eigenvalues of Htoto are then the same 
as ()126p . But the bases of unperturbed SESR can be taken as 

Ne 

\ip"^nE) = \ij"^)(E)<^\nk). (127) 

k=l 

Case three: The Hamiltonian split is 

-fftoto = or constant, Htoti ^ fJ-c^i Ie + cr^ <^ (B^e + B^e) ■ (128) 

The bases of a selected SESR are just the natural bases I^s^-e) defined in (lllip . Then, we use our Hamiltonian 
redivision skill to obtain 

Hito = i/zurck, i?t'otl = M<^S + cr^ B^E- (129) 

The corresponding eigenvalues of H^^^q is given by (|114p . 

Case four. The Hamiltonian split is the same as ()128p . But the bases of the unperturbed SESR can be chosen as 



l^sXnJ = l"s)®(g)lx"'=)- (130) 

k=l 

Then, we use our Hamiltonian redivision skill to obtain 

H^to = iJzurek + ctI(E) B^e, H'^^^ = A't^s ■ (131) 

The corresponding eigenvalue is 

Ne 

E'n.nE = Ksnrn.- = ("1)"' ^.^-IT'-Y,. (132) 

fc=l 
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It must be emphasized that the four kinds of choices on the SESRS aim at the different preconditions if the cut-off 
approximation of perturbation is necessary. Cases one and two are used to the preconditions that fi ^ Zi^ and/or 
/i ^ Xk; that is, the transverse field fi is strong. Case three is chosen when Zk ^ IJ. and Zf. ^ Xk- In other words, 
two transverse fields are weak. Case four is suitable to solve the problem under 3> /i and/or X^ 3> This means 
that the transverse field fi is weak. 

It is easy to see that in cases one and two there are two degenerate subspaces with Ne dimensions, which cannot 
be completely removed via the usual diagonalization procedure of the degenerate subspaces and our Hamiltonian 
redivision. However, the conditions that degeneracies happen are Smsn^- For case one 



msrriE.nsriE _ r 
Hi ~ "ms(l-ns) 



while for case two 



Oms(l — ns) 



E 

k=l 



i-iy 




Lk=l 



Ne 

E 

fc=i 



Ne 



1 = 1 



11'^" 



XkS, 



'mfc(l-nfc) 




(133) 



(134) 



This implies such a fact that in both case one and case two t?]^ 



levels, that is 



msmE,nsnE 



are vanishing between any two degenerate 



^msiis — 0. Therefore, our improved scheme of perturbation theory can work well. However, 



note that the preconditions that ii ^ Zk and/or /i '3> Xk in cases one and two are the same, we prefer to use the 
choice of case one because its calculation is easier than case two in our improved scheme of perturbation theory. 

As to case three, we also can not completely remove the degeneracies via the usual diagonalization procedure of 
the degenerate subspaces and our Hamiltonian redivision. The conditions that degeneracies happen are solutions of 
the following equation 



AfE 

E^4(-ir 

k=l 



while the off-diagonal elements of the perturbing Hamiltonian matrix are 



msmE,ns"E _ ,,s; 



Ne 

ms(l-"s) W ^" 
k=l 




™fc(l-"/c) 




(135) 



(136) 



It is clear that we can not guarantee, in general, that g™smE,nsnE vanishing between any two degenerate levels. If 
the result is indeed so. This SESR is not a good choice because the remained degeneracies will result in the difficulty 
to use the usual perturbation theory and complication in our improved ones. If we still intend to use the cut-off 
approximation of perturbation, the results from case three will not be satisfied enough if the evolution time is long 
enough, because Ai (t) or An (t) has the extra terms proportional to the evolution time that can not be simply absorbed 
to the exponential power for / > 2 in our views. 

Fortunately, that case four can be covered the precondition that Zk ^ /i and Zk S> Xk in case three. Hence, we 
give up the choice of case three and only use the SESR in case four. Actually, the above problems originally motivate 
us to consider how to choose the appropriate SESR for open systems, which has been seen in Sec. [Hi 

It is easy to get that the conditions that degeneracies happen in case four are solutions of the following equation 



Ne 



J2 Yk [(-1)"='+™'= - (-l)"s+n^] = 0, 



k=l 



while the off-diagonal elements of the perturbing Hamiltonian matrix are 



Ne 



msmE,ns«E 
i/l 



k=l 



Our interesting task is to seek for the conditions that degeneracies happen when g" 
iTik = rik for any k, into Eq. (|137p and rewrite it as 



(137) 
(138) 

msmE,«s«E Hence, we substitute 



Ne 



E^^-(-i) 



.fc=i 



[(-1)™^ - (-1)"^] = 0. 



(139) 
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Its solution is ms = ns unless the exception X^fcSi '^ki—^)"^'' = 0. However, this exception is not valid limiting our 
problem to open systems because it means that Hj-oto — from Eq. (|132p . or equivalently, the total Htot becomes 
/X(Ts d) Ie- Again jointly considering it and the expression (|138p of g^sniE .nsuE ^^^^ four, we obtain the conclusion 
that g^smE,nsnE ^^^^ indeed vanishing between any two degenerate levels. 

In the following discussion, we only focus on case one with the strong transverse field /i and case four with weak 
transverse field /i in order to illustrate our exact solution and improved form of perturbed solution more simply and 
better. 

Let us define 



Ne 



(140) 



ni,n2,-"— nE— 

Thus, for case one, we can rewrite the off-diagonal elements of the perturbing Hamiltonian matrix as 

711 

9i 



— Oms(l-ns)/™E""iE"B- 



Substituting it into the definition of GngnE ; obtain 



G 



(2) ^ JriE ^(3) ^ Q r;(i) 

2^ ' ^ns"E ' ^nsn 



Hence, we have 



Er, 



(-l)«^/i 



1 



1 fL 



2 \ 2fi 

Similarly, for case four, from Eqs. (|132p and (|138p it follows that 



1 (fuE 



G 



^(-1)-^ ^ ^_(-irv 



2/, 



8f3 



Hence, we have 



(141) 



(142) 



(143) 



(144) 



(145) 



En 









( ^ \ 











(146) 



It is clear that there is a corresponding relation between the case of the strong transverse field /i and the case 
of weak transverse field /i, that is, their perturbed solutions will be the same under the exchanging transformation 
/i ]n-E- Hence, we only write down, for case one, the zeroth, first and second order parts of total system density 
matrix at time respectively 



mg,mE— ns,nE— 



i B, 



.s™E-5-.s"E)v™f'"-"^"-(o)iV'r)(V'r I ® ix'"-)(x"-i, 



(147) 



ms,mE— ns,nE— 



msmE.rasnE 
^Ptot 



E E 

(0)IC)(V's""1®lx"^>(x 

E E fe-'(^"s-E-^"s"E)*_ 



2fj. 



mQ,m-E—0 ng ,nE— 



2/x 



^Ptot 



(148) 
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ms,mE— ns,nE— ^ 



/tie 

2^ 



/toe 

2^ 



E E (^-' 



ms,mE— ns,nE— 



') 



(i-„s)„ 



"E msmE,nsnE 



(0)|V'^™^)(V's""1«lx'"^)(x"^|. 



(149) 



It is easy to give the solution of the reduced density matrix of the open system up to the improved form of the second 
order approximation by tracing off the environment space, that is 



(150) 



For a given initial state, this trace is very easy to calculate and the explicit form of of the open system solution is 
obtained. Then we can discuss the decoherence and entanglement dynamics according to the methods in, for example, 
[igL [13] , and they are arranged in our forthcoming manuscript (in preparing) [l8l | . 



IX. DISCUSSION AND CONCLUSION 



This paper studies open systems dynamics, which is the third in our serial studies on quantum mechanics in general 
quantum systems. Its conclusions are obtained based on our previous two works (l^. [l3j. 

It must be emphasized that we study open system dynamics according to the "the first principle" , that is, the 
Schrodinger equation and the von Neumann equation, and we do not consider the phenomenological methods and 
theories. For generality in theory, we obtain the exact solution of the open system without using any approximation. 
The deduction of our exact master equation only uses the factorizing initial condition. Particularly, we derive out 
our perturbed master equation and its improved form, but we give up all of approximations used in the traditional 
methods and formulism except for the factorizing initial condition. It is very interesting that we get the Redfield 
master equation without using the Born-Markov approximation. This implies the Born-Markov approximation is 
unnecessary based on our results. 

A simple but key idea to obtain our exact solution of open systems is an appropriate choice of the SESR. In 
fact, it closely connects with the Hamiltonian redivision skill [l3|. In Sec. IVIIH we have clearly stated its reasons. 
Originally, the aim that we propose this idea is to break the accustomed choice of + H^, and build a picture to 
allow the interaction between the open system and its environment into our unperturbed representation. This makes 
the Hamiltonian redivision skill to look like more natural. 

Our exact solution and master equation of open systems are general and explicit in form because all order approx- 
imations of the perturbing Hamiltonian not only are completely included but also are clearly expressed, although it 
is an infinite series. In special, they are in c-number function forms rather than operator forms. This means that 
they can inherit the same advantage as the Feynman path integral expression. Moreover, they are power series of the 
perturbing Hamiltonian like as the Dyson series in the interaction picture. This implies that the cut-off approximation 
of perturbation can be made for the needed precision of the problems. 

Based on our improved scheme of perturbation theory, the improved forms of perturbed solution and perturbed 
master equation can absorb the partial contributions from the high order even all order approximations of perturbation. 
Therefore, we can say that our open system dynamics is actually calculable, operationally efficient, conclusively more 
accurate. 

In order to extend our method, we also discuss Milburn model of open systems. In fact, from our point of view, 
Miburn model of dynamics should be applied to so-called Milburn-type closed quantum systems made up of the 
interesting open system and its near environment. A Milburn-type closed quantum system is not really closed system 
from the view that a really closed system must obey the von Neumann equation. If one cannot know the Hamiltonian 
of its remote environment and the interaction form between the interested system and its remote environment, Milburn 
model of dynamics might be a choice scheme to study this kind of open systems. In the above sense, the extra term 
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with ^0 multiplier in the Milburn equation represents the influence from the remote environment. Obviously, when 
^0 = 0, Milburn dynamics back to von Neumann dynamics. This implies that the (very) remote environment can be 
ignored. We obtain the exact solution that can provide a general tool to investigate those interesting and complicated 
open systems when the environment model is partially known. However, there a free parameter 6q in the Milburn 
model. It is still not been given by the theory, but it should be able be known by the experiment if we think the extra 
term in the Milburn dynamics is resulted in by the remote environment. 

Note that our open system dynamics is derived from the first principle, our open system dynamics is not applicable 
to the cases that ones do not clearly know the Hamiltonians of the open system, its environment and the interaction 
between the system and the environment unless at this time Milburn model is suitable. How to relate with some 
phenomenological theory of open systems will be done in the near future. 

As examples, Zurek model of two-state open system and its extension with two transverse fields are studied, 
respectively, in the strong and weak fields acting on the system. We specially display how to choose the appropriate 
SESR. They indicate that our open system dynamics is a powerful theory and tool. We are sure that our open system 
dynamics can be used to more open systems since its generality and clearness, and the calculations are simpler and 
more efficient, the results are more accurate and more reliable than the existed scheme. 

In summary, our results can be thought of as theoretical developments of open system dynamics, and they are 
helpful for understanding the theory of quantum mechanics and providing some powerful tools for the calculation 
of decoherence, entanglement dynamics, quantum dissipation, quantum transport in general quantum systems and 
so on. Together with our exact solution and perturbation theory [T^ . [Tsj . they can finally form the foundation of 
theoretical formulism of quantum mechanics in general quantum systems. Further study on quantum mechanics of 
general quantum systems is on progressing. 
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